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Direct numerical simulations (DNS) of particle-laden turbulent flow in straight, mildly 
curved and strongly bent pipes are performed in which the solid phase is modelled as 
small heavy spherical particles. A total of seven populations of dilute particles with dif¬ 
ferent Stokes numbers, one-way coupled with their carrier phase, are simulated. The 
objective is to examine the effect of the curvature on micro-particle transport and ac¬ 
cumulation. It is shown that even a slight non-zero curvature in the flow configuration 
strongly impact the particle concentration map such that the concentration of inertial 
particles with bulk Stokes number 0.45 (based on bulk velocity and pipe radius) at the 
inner-bend wall of mildly curved pipe becomes 12.8 times larger than that in the vis¬ 
cous sublayer of the straight pipe. Near-wall helicoidal particle streaks are observed in 
the curved configurations with their inclination varying with the strength of the sec¬ 
ondary motion of the carrier phase. A reflection layer, as previously observed in particle 
laden turbulent S-shaped channels, is also apparent in the strongly curved pipe with 
heavy particles. In addition, depending on the curvature, the central regions of the mean 
Dean vortices appear to be completely depleted of particles, as observed also in the par¬ 
tially re-laminarised region at the inner bend. The turbophoretic drift of the particles 
is shown to be affected by weak and strong secondary motions of the carrier phase and 
geometry-induced centrifugal forces. The first and second-order moments of the velocity 
and acceleration of the particulate phase in the same configurations are addressed in a 
companion paper by the same authors. The current data-set will be useful for modelling 
particles advected in wall-bounded turbulent flows where the effects of the curvature are 
not negligible. 


1. Introduction 

Most of the flows observed in nature and industrial applications are turbulent, and 
are often carrying dispersed particulate phases. Some typical examples are the atmo¬ 
spheric boundary layer transporting volcanic ash and pollen particles, or the flow in 
engine combustion chambers with suspended droplets. In the context of turbulent flows, 
the dispersion process of heavier-than-fluid particles predominantly results in two pecu¬ 
liar phenomena: small-scale clustering and turbophoresis. Small-scale clustering affects 
the particle distribution resulting in the loss of spatial homogeneity since the parti¬ 
cles aggregate in specific regions characterised by higher values of the turbulent kinetic 
energy dissipation rate; this phenomenon appears in all kinds of turbulent flow laden 
with heavy particles. On the contrary, turbophoresis is typical of turbulent wall-bounded 
flows and occurs as a pronounced particle accumulation in the near-wall region. These 


f Email address for correspondence: azad@mech.ktli.se 


2 


A. Noorani, G. Bardina, L. Brandt, P. Schlatter 


two phenomena are not completely independ ent but it has been sh own that small-scale 
clu stering' induces turbophoresis in w all flows (Sardina et a/.l2012Q, l. Recent publications 
by Balachandar fc EatonI ( 2010l l and Toschi &: Bodenschatz ( 2009 1 provide a review of 
the research activities in the area of dispersed phase in turbulent flows. 

The origin of the preferential particle concentration resulting in small-scale cluster¬ 
ing is still being studied. It is usually explained in terms of the combination of parti¬ 
cle inertia and turbulent vortical structures: Inertia prevents particles from following 
the fl uid trajectories leading to preferential conce ntration outside of the vortical re¬ 
gions ( Squires fc Eatonlll99ll : IWang fc Maxev^lI993^ . Other explanations have also been 
proposed for small-scale clustering since this peculiar phenomenon is also observed in 
random fl ows without the co herent vortical structures that usually appear in turbu- 
lence (see Mehli g et a/.l [2005). Among these, the sweep- stick mechanism proposed by 


Goto fc Vassilicos ( 2008h and Coleman fc Vassilicod ( 2009l l suggests that the inertial par¬ 


ticles tend to preferentially accumulate in regions characterised by a stagnation point of 
the flow acceleration field. 

On the other hand, turbophoresis is essentially an inertial p article drift towards the 


walls driven by turbul ent fluctuatio ns (ICaporaloni et al. 1197511 . This process has been 


theoretically studied by ReeksI ( I983ll : Young fc LeeminJ i I997tl . and sever al experiments 


have b ee n conducted to observ e and quantify this near-wall accumulation (|Kaftori et al. 


iggSol ffc: Ninto fc GarciallT996ll . Today, the most common way to address the problem is 


by means of direct numerical simulation (DNS ) in classical channel or pip e-flow configu¬ 
rations at moderate Reynolds number (see e.o. ISoldati fc Marchioli ( 2009h for a review). 
In particular, it is assumed that turbophoresis is strongly linked to the vortical structures 
of near-wall turbulen ce. Sweep and ejection e v ents affect the particl e trans fer towards the 
wall as described by Rouson fc EatonI ( 2001 1. Marchioli fc Soldati ( 2002h observed that 
the inertial particles tend to preferentia lly oversample t he low streamwise velocity regions 
in the buffer layer, i.e. ejection events. ISardina et all (12012)71 1 further showed a link be¬ 
tween this preferential particle localisation in the l ow speed streaks a nd the steady-state 
concentration profiles of the particles. According to Milici et a,l. 1 2014h turbophoresis may 
be completely suppressed when the streaky patterns of wall-turbulence are destroyed e.g. 
by adding wall roughness to a channel flow. The turbophoretic motion of the particles is 
also present in boundary-layer flows since the turbulent struct ures close to the wall are 
essentially simi lar to those of the canonical channel flow (see Sardina et al. 2012lJ l. A 
recent study bv ISikovs^ ( 20I4h investigates the stochastic behaviour of inertial particles 
near the wall by means of matched asymptotic expansions. This study suggests that the 
particle concentration in the viscous sublayer of wall-bounded turbulent flows exhibits a 
power-law singularity and the corresponding exponent depends on the particle relaxation 
time. 

Internal flows relevant to industrial applications often contain parts with geometrical 
complexities such as curved sections or converging/diverging compartments, which in 
general produce secondary flows, separation and other complex flow phenomena. Despite 
some previous efforts, particle motion has mostly been studied in the context of canonical 
wall-bounded flows. Given the strong link between turbophoresis and turbulent near-wall 
structures, it is fundamental to investigate the particulate phase behaviour in complex 
geometries where a slight modification of the geometrical characteristics of the flow can 
induce a significant change in the vortical wall-dominated structures and consequently 
in the turbophoretic drift and particle clustering. The current research, therefore, will 
examine the effects of curvature on particulate dispersion in bent pipes in comparison to 
the canonical (and more studied) straight pipe. More specifically, turbulent flow in curved 
pipes is occurring in many typical engineering applications such as heat exchangers. 
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chemical reactors and pipeline systems. In these typical heat and mass transfer systems 
the flow is commonly laden with a dispersed phase e.g. membrane ultrafiltration hydraulic 
systems or mixing devices applied in pharmaceutical industries. 

The imbalance between the cross-stream pressure gradient and geometry induced cen¬ 
trifugal forces causes a secondary motion orthogonal to the main flow direction, which is 
characterised by the formation of a pair of counter-rotating vortices, the so-called Dean 
vortices. These vortices force the streamwise velocity to be distributed non-uniformly in 
the cross section of the pipe such that its maximum is deflected towards the outer wall. 
Secondary motion that is geometry or skew-induced appears both in laminar and turbu¬ 
lent flows unlike the stress-induced secondary flow that is solely formed as a result of the 
local variation of the Reynolds stresses in turbulent regimes. The former type is usually 
referred to as Prandtl’s secondary flow of first kind and the latter is called Prandtl’s 
secondary flow of second kind which can be observe d, for instance, in the turbulent flow 
in straight ducts with non-circular cross-section (see Bradsha^ 1987ll . 

In general, curved pipe geometries can be classified as: (a) Helically coiled tubes where 
the flow is fully developed inside the curved geometry; (b) Spatially developing bends 
such as {/-bends or elbows (90° bends) where a fully developed flow from a straight 
pipe enters the bend. Generally the pitch of the resulting coil in category (a) causes 
an additional force which is the torsion acting alongside the centrifugal force on the 
fluid flow. However, when the coil pitch is small compared to the coil diameter, which 
is the case in most practical applications, the influence of the torsion i s negligible wit h 
respect to the effect of the curvature {see iManlaoaz fc Churchilll 1981 : Germanolll982 b 
This reduces the geometry to that of a toroidal shape. The latter category (b) and in 
particular bends with small curvature where the entry flow region is short compared to 
the length of the pipe in the curved section can also be modelled as toroidal pipes. This 
abstraction - infinitely long bent pipe configuration - thus provides a unique opportunity 
to isolate the effect of the curvature on the turbulent characteristics and also to study the 
influence of the centrifugal forces and the secondary motion on the near-wall features. 

Due to the complexity of the flow field in curved pipe configurations very limited efforts 
have been made i n the past to unfold the intrinsi c dynamics related to the turbulent flow 
in toroidal pipes. Boersma fc Nieuwstadt ( 1996ll employed large-eddy simulation (LES) 
to examine the effect of curvature on the mean flow and fluctuations. iHiittl fc Friedrich 

S performed the first DNS in the framework of Germano’s coordinates ( Germano 
to test the curvature and torsion effect on turbulence at a friction Reynolds number 
Rcr = 230, b ased on the az i muth ally averaged friction velocity and pipe section radii. 
Only recentlv. lNoorani et all ( 2013ll performed DNSs with high enough Reynolds numbers 
for the flow to remain turbulent even in strongly curved pipes. These authors computed 
the full Reynolds stress budget at various curvature configurations for the first time 
and have shown that the turbulent flow in the inner side of the strongly curved pipe is 
highly damped while the outer bend remained fully turbulent with essentially unchanged 
near-wall dynamics. 

A dispersed phase transported in turbulent curved pipes is subjected to the wall- 
dominated turbulence of the carrier phase, its secondary motion due to the Dean vor¬ 
tices and volumetric geometry-induced centrifugal forces. These effects may have a large 
impact on the turbophoretic motion of the particles. Investigations of a dispersed phase 
in turbulent flow in toroidal pipes are difficult to per form. Therefore, num erical or ex¬ 
perimental data are very rare in the literature (see Vashisth et all 2008h . Examining 
membrane filtration modules, as an example, shows the reduction of near-wall concen- 
tration of the particles in curved ge ometries compared t o the straight configuration (see 
Winzeler fc Belfortlll993l l. Recently. IWu fc Young ( 2012 ) conducted a set of experimental 
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measurements and theoretical calculations considering the inertial particle dispersion rate 
in spatially developing mildly curved ducts (with rectangular cross-section). In compari¬ 
son to the straight configuration, they reported a dramatic increase in the accumulation 
rate of particles at the outer bend surface. Their study showed that large heavy iner¬ 
tial particles are mostly driven by centrifugal forces rather than turbophoresis even in 
a mildly curved configuration. From a numerical point of view, some attempts to study 
particle transport i n 90° pipe bends have be en recently conducted by means of LES 
( Breuer et a^. boO^ iBerrouk &: Laurencell2008[l . Applying an artifi cial swirl in a spatially 
developing turbulent straight pipe simulation Zonta etall ( 2013h studied the effect of 
centrifugal forces along with turbophoretic drift on near-wall particle accumulation. 

In anot her study to exam i ne th e part icles behaviour in the presence of geometrical 
curvature, iHuang fc Durbiii ( 2010l 12012 1 performed DNS of particulate dispersion in 
a strongly curved S-shaped channel. They examined the erosive effect of wall-particle 
collisions at the channel walls and reported that the maximum erosion rate is dramatically 
increased for heavier particles. They also observed a plume of particles bouncing back 
from the outer wall of the curved section and building up a hi gh-particle concen t ration 


. _ iig£ _ _ 

region adjacent to the wall, which is named as ’reflection Zar/er I IHuang fc DurbinI ( 2012 1 
further presented a simple model to explain the oscillatory behaviour of the wall collision 
for heavy inertial particles. 

The current study is devoted to the complex dynamics of particle dispersion in the 
presence of the curvature where a stably secondary motion and centrifugal forces acting 
on the particles co-exist. Fully developed turbulent pipe flows at three different curva¬ 
tures laden with inertial small heavy particles are simulated. In these configurations, 
the transport and accumulation of these dilute micro-size particles are investigated and 
addressed in the current paper. The first and second-order moments of the velocity and 
accelerat ion of the parti c ulate phase in turbulent bent pipes are addressed in a companion 
paper bv iNoorani et al\ ( 2015 1. In the pages that follow, the computational methodology 
for simulating the carrier and dispersed phase is explained in along with a discussion 
of the basic flow features in 1 13.11 The instantaneous particle distribution in the pipe, 
particle trajectories, and Eulerian statistics of the Lagrangian phase, in particular con¬ 
centration maps, are discussed in 93.21 Conclusions and an outlook are given in the final 
section El 


2. Computational Methodology 

2.1. Flow configuration and governing equations 

In order to provide a clear overview of the computational domain employed for the DNSs 
a schematic of a curved pipe is given in figure [TJ The pipe consists of a part of a torus 
whose centre is located at the origin of the Cartesian system {X,Y,Z). Relative to this 
system the toroidal coordinates {R, s, C) of the pipe together with local (in-plane) poloidal 
coordinates {r,0) are shown in figure [T] (a). The equatorial mid-plane of the pipe which 
is also the plane of symmetry is indicated as vertical cut of the cross-section of the pipe 
at 0 = 7r/2. The horizontal cut of the pipe section at 0 = 0 is displayed in figure [T] (b). Of 
particular relevance in this configuration is the curvature parameter k that can be defined 
as Ra/Rc where Ra is the radius of the pipe cross-section and Rc is the major radius of 
the torus at the pipe centreline (shown in figure [T]c). Generally k distinguishes between 
mildly (weakly) curved pipes (/t 0.01) and strongly curved configurations (/t > 0.1). 

All simulations are performed in Cartesian coordinates where the incompressible Navier- 
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Figure 1. (a) Side view of curved pipe geometry and its associated Cartesian coordinates, 
{X, Y, Z), with embedded toroidal coordinate system, {R, s, (^). The in-plane polar coordinates, 
(r, 9), relative to the cross-section Cartesian triads, {x, y), of the pipe-section are also illustrated 
at the left panel, (b) Cut-away view of a curved pipe, with locations of the horizontal and vertical 
cross-sections displayed at 9 — 0 and 0 = 7r/2 respectively, (c) The schematic front view of the 
curved configuration with toroidal and in-plane poloidal coordinates. 


Stokes equations for the carrier phase read in non-dimensional form 


V • M = 0, 

— -f (m • V)m = - Vp -h . 

ot neb 


( 2 . 1 ) 

( 2 . 2 ) 


Here, u indicates the velocity vector of the fluid flow, and p is the pressure. Rcb stands 
for the bulk Reynolds number dehned as 2RaUb/v', Ub is the bulk velocity, and u is the 
kinematic viscosity. The Dean number is often used in the literature to characterise the 
secondary motion and is defined as Dcb = Reby/l^ for the current flow configurations 


( Berger et al.lll983 1. 


2.2. Numerical approach for the carrier phase 
The present DNSs are performed using the spectral-e lement code nekhOOO. This mas¬ 
sively parallel code is developed by Fischer et a,l. ( 2008h at Argonne National Laboratory 
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Figure 2. (a) A quarter of a cross-sectional plane pertaining to the pipe mesh, (b) Grid represen¬ 
tation in the equatorial mid-plane of the toroidal pipe. Here, element boundaries and Gauss-Lo- 
batto-Legendre (GLL) points are shown in the darker region. 


(ANL). The physical domain is divided into a number of hexahedral local elements on 
which the incompressible Navier-Stokes equations are solved by means of local approx¬ 
imations based on high-order orthogonal polynomial basis on Gauss-Lobatto-Legendre 
(GLL) nodes. This provides spectral accuracy with geometrical flexibility applicable to 
problems with moderately complex geometries. The spatial discretisation is obtained by 
means of the Galerkin projection, applying the Pat — Pw sche me where the velocity an d 
pressure spaces are represented by the same polynomial order ( Tomboulides et '^ Il997ll . 
The temporal discretisation uses a 3’’'^ order mixed Backward Difference/Extrapolation 
(BDF3/EXT3) semi-implicit scheme. The walls are set to no-slip boundary condition 
while periodicity is applied for the end sections of the pipe domains in the axial direction 
as the flow is assumed to be homogenous in this direction. The initial condition for all the 
simulations are set to be the laminar Poiseuille profile of the straight pipe. This initial 
flow is perturbed with low-amplitude pseudo-random noise which is then evolved in time 
until the turbulent state is statistically fully developed in the pipe; this corresponds to 
approximately 200 integral (convective) time units from the start. The nondimensional 
convective time unit is defined as tub/Ra- The driving force of the flow is adapted in such 
a way to keep the mass flow rate fixed throughout the simulations. This is performed 
implicitly in the time-integration scheme. 

Note that the nu merical setup for the carrier phase is similar to the one used by 


Noorani et al\ (1201311 for their highest Reynolds number; in particular the mesh topology 


and resolution are the same. Nevertheless, we provide a short overview of the main char¬ 
acteristics here for completeness. Figure [2] (a) displays a quarter of the cross-sectional 
grid that is generated for the pipe section. The two dimensional mesh is uniformly ex¬ 
truded in the axial direction to generate the hexahedral elements of the 3D straight pipe 
mesh. The required grid for curved/toroidal configurations is obtained using an analytical 
morphing of the straight pipe grid. For all the cases, 8^ GLL points (polynomial order 
7) resolve the grid inside each individual element. The resulting mesh in the equatorial 
section of the bent pipe is sketched in figure [2] (b). 

It is known that the critical Reynolds number for transition to turbulence is generally 
higher in curved pipes compared to the straight configuration. The present simulations 
are carried out at a fixed Rcb = 11700 such that the flow in all curved configurations is 
above the critical value given by Ito ( 1959ll . The doma i n leng th of all configurations is 
fixed to be 25i?Q along the pipe centreline. iGhin et al\ ( 2010!) studied the effect of the 


streamwise periodic length on the convergence of turbulence statistics in the straight 
pipe. Their study concludes that a pipe length of 87ri?a is sufficient for the flow statistics 
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Rei, ^ of elements # Grid points Ar'*' A{Ra6)'^ As^ 


11700 237120 121.4 x 10® (0.16,4.70) (1.49,4.93) (3.03,9.91) 

Table 1. Resolution details for the present pipe flow meshes. All quantities are based on the 

straight pipe configuration. 


Reb,D K Rct Averaging time frame (tUb/Ra) 


11700 0.00 360 3000 

11700 0.01 368 3000 

11700 0.10 410 3000 

Table 2. Simulation parameters of the present study. The time t is normalised by Rajub- The 
non-dimensional friction Reynolds number Rct is defined as uPRajv. 


to be independent of the periodic domain length; our choice of the domain length is 
based on these findings. In the straight configuration the grid spacing is designed to 
have 14 grid points placed below y’*' = (1 — rlRfly^ = 10, measured from the wall, and 
A(i?a0)max ^ ^ and < 10. These are measured in wall units where the viscous 

length scale is based on the straight pipe friction Reynolds number (the superscript + 
indicates the viscous scaling). The resolution of the grid for the curved pipes is chosen 
to be the same as that of the straight configuration. The characteristics of the grid are 
summarised in table [T] 

After the fluid flow reached the statistically stationary state, statistics of the turbulent 
flow are computed based on averages (denoted by (•)) in time (t), and along the axial 
direction (s) which is the homogeneous direction of the flow. The statistical sampling 
is started approximately after 200i?a/itb to make sure the flow is fully developed and 
settled in the turbulent state. The averaging time frame and other relevant simulation 
parameters are listed in tabled Note that for all cases 3000 convective time units were 
included in the temporal average, which proved to be necessary for obtaining converged 
statistics for the Lagrangian phase; for the flow statistics a shorter time frame would have 
been sufficient. Further details on how to obtain statistical averages via tensor transfor¬ 
mations unde£j;otatiOT from Cartesian coordinates to toroidal/poloidal coordinates are 


given in Noorani et all 1 2013h . The verification of the method and the validations of the 


simulation setup alon g with statistical analysis for the carrier phase is also performed by 


Noorani et al. ( 20131) . 


In the curved configurations the mean friction velocity Ur is defined as \jTw,totlpp 

where rTTflot indicates the total mean shear stress over the wall and p/ is the fluid density. 

- 2 - 2 

The mean shear stress can be obtained as ^{{tui,s) + (Ru.e) ) where overbar denotes 
averaging along the circumference of the pipe section of the streamwise mean shear 
stress {tw,s) and its azimuthal component {twp)- Note that for the straight pipe, due to 
symmetry {Tw,e) vanishes. The resulting is used to compute the viscous scaling. 


2.3. Numerical approach for the particulate phase 

In order to characterise the dispersed phase one commonly defines the volume fraction of 
particles with respect to the fluid as = NpVp/Y where Np is the number of particles. 
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Yp is the volume of one particle and V is the global volume occupied by the carrier 
phase. The density ratio is defined as Pp/pf, Pp being the particle density. According to 
Elghobasliil ( 199ll . [l994 1 these two parameters ($„ and Pp/pf) determine the level of in¬ 
teraction between the dispersed and carrier phase. Considering sufficiently small volume 
fraction < 10“® and large enough density ratio (around ^ 10^) the dynamics of the 
particles is essentially governed by the turbulent carrier phase. In fact, for dilute small 
particles the feedback for ce acting on the fluid phase and the inter-particle collisions can 
be considered negligible ( Balachandar &: Eatonll201Cll) . In this case the inerti al particles 
are one-way coupled with the carrier (Eulerian) field. Maxev Sz Rilev ( I983ll rigorously 
derived the equation of motion for a spherical rigid particle immersed in a non-uniform 
flow. The resulting equation of motion includes various forces acting on the individual 
particles. These derivations have be en subjected to a n umber of corrections and some¬ 
times additional terms are included ( Crowe et Qi.ll20IlI) . Among the different forces the 
only significant contributions to driv e a heavy point particle are th e aerodynamic Stokes 
drag and the gravitational force Isee lElghobashi fc Truesdeilll992l) . 

The current research is concerned with micron-sized dilute heavier-than-fluid particles 
that are assumed to be smaller than the smallest spatial scale of the flow. With these 
assumptions, the particle equation of motion reads 


dvp 

dt 

dxp 

dt 


(I + 0.15i?e°-687) , 

btb 




(2.3) 

(2.4) 


where Xp is the particle position, Vp and u{xp,t) are the particle velocity and the fluid 
velocity at the particle position, respectively. The non-dimensional bulk Stokes number 
Stb is defined as the ratio of particle relaxation time Tp to bulk flow residence time 
{Ra/ub), and characterises the particle response to the fluid. It can be computed as 
TpUb/Ra where Tp = ppd^/{18pfv)-, dp is the particle diameter. The Stokes number as 
dimensionless form of the particle relaxation time can also be expressed considering the 
viscous (turbulence) time units such that St~^ = TpiGB/v. The Stokes number based 
on the friction velocity {St'^) is more meaningful in the straight pipe/channel where 
near-wall turbulence is responsible for nonuniform inertial particles dispersion. 

Equation (I2.3|l essentially expresses the particle acceleration (op) due to a steady- 
state drag force acting on the spherical solid particle in a uniform velocity field. The last 
term on the rig ht-hand side is the non-lin ear correction due to particle finite Reynolds 
number effects ( Schiller &: Naumannlll93^ . The particle Reynolds number Rcp is defined 
as \vp—u{xp, t)\dp/v. For the current investigation, the effect of gravitational acceleration 
will be neglected to avoid losing generality and to be able to isolate the effect of the 
centrifugal force and the secondary motion of the Dean vortices on the turbophoretic 
behaviour of particles. 

Unlike the carrier phase which is simulated on a set of fixed grid points (Eulerian 
framework), the particles are tracked freely inside the computational domain (Lagrangian 
framework). The accuracy of the interpolation method is significant for the correct eval¬ 
uation of the hydrodynamic forces and for inst ance it has direct conseq uences for the 
acceleration spectrum of the fluid particles ((see Ivan Hinsberg et 0,1.1120131 1. A spectrally 
accurate interpolation scheme, hence, is applied in the current study to evaluate fluid ve¬ 
locity at the centre of the spherical particle. Since in each time step the particle tracking 
is performed after updating the fluid phase, the Legendre basis functions and the corre¬ 
sponding expansion coefficients are already available for each element. Therefore, without 
losing the accuracy, the fluid properties (velocity, pressure, etc. ) can be determined in 
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arbitrary positions inside each element. The order of accuracy of this interpolation, thus, 
is equal to the order of spectral element method (7 in the current simulations). The time 
stepper implemented for the particle tracking is the classical 3’’'^ order multistep Adams- 
Bashforth {AB3) scheme. The time step of the carrier phase solver (At) is fixed such that 
the convective CFL condition is always below 0.5. The time step of the particle solver is 
also set to be equal to the flow solver time step. As for the velocity field, re-starts of the 
simulations are based on storing all necessary previous times in order to avoid reducing 
the temporal order. 

The particle interaction with solid walls is treated as elastic collisions. A restitution 
coefficient can be defined as e = \vn, 2 lvn,i\] Vn ,2 and Vn,i are the wall-normal pre- and 
post-collisional particle velocities. In a purely elastic collision (as used in the current 
study), e equals unity and the total kinetic energy is conserved in the collision process. 
It has been shown previously that inelastic reflection from the wall has a n egligibly 


small effect on the dynamics o f small inertial pa r ticles in turbulent flows (see iLi et al 


l200lh . However, the analysis of iHuang fc Durbin ( 20I2h suggested that in the presence 
of centrifugal forces the amount of collisions and their significance will increase especially 
for particles with larger St. Applying an inelastic collision scheme might therefore slightly 
alter the turbulent particle structures. Similar to the carrier phase, periodic boundary 
conditions are used for the dispersed phase in the axial direction. 

The implemented particle tracking scheme is validated for a classical particle-laden 
turbulent channel flow. For that purpose, DNS at Rer = 180 (based on the half-height 
of the channel) in a stan dard domain ( 47r x 2 x 4/37r) was performed and compared 
with the existing data of Sardina et al.l (|2 0I2(Jl. The latter s imulation was performed 
using the fully spectral code SIMSON ( Chevalier et al ] [200i . A total of four particle 
populations is considered with Stokes number (S't+ = 0,1,10,100) where Pp/Pf is set to 
770. The two simulations share similar set-ups in terms of initial conditions, start and 
end time of the particle tracking and statistical sampling time frames. The normalised 
particle concentrations in the wall-normal direction, shown in figure [H show an excellent 
agreement between the data obtained with the two different codes. 


3. Results and analysis 

In this section, results from the simulations detailed in lj2]are discussed. In a first step 
in EH the carrier phase is discussed with some detail, however the main focus is on 
characterising the particle phase in (13.21 


3.1. Carrier phase 

It is crucial to have a clear overview of the carrier phase turbulence before analysing the 
complex dynamics of a particulate phase in bent pipes. To this end the main features 
of the turbulent flow in these configurations are initially described and discussed in this 
section. A mo r e det ailed analysis of the turbulent curved pipe flows can be found in 


Noorani et all ( 20131) which uses the same numerical setup. 


The qualitative modifications of the flow when increasing the curvature parameter are 
illustrated in figure |4] where instantaneous cross-sectional views of the axial velocity are 
shown in the equatorial mid-plane and at an arbitrary axial position. Evidently, as the 
curvature parameter increases the bulk flow is deflected further towards the outer bend, 
while the maximum of the axial velocity remains in the centreline of the straight pipe. 
This leads to a reduction in turbulence activity with increasing curvature. In the strongly 
curved pipe, this turbulence depression is large enough to almost re-laminarise the flow 
in the inner bend (c./. figure |6]c). 
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Figure 3. Normalised statistically steady-state particle concentration distributed in the wal¬ 
l-normal direction of the tur bulent channel . . = 0,-= 1, - St'^ = 10, 

-= 100 ; □ d ata by ISardina et al\ l|2012af ) with the pseudo-spectral code SIMSON 

llChevalier et fflZ.ll2007fl : • present nekSOOO simulation. 


With increasing curvature the centrifugal forces increase. From figure [5] it is notice¬ 
able that the maximum of the mean streamwise velocity deflects more towards the outer 
bend with increasing k. The magnitude of the averaged in-plane velocity + (ug)"^) 

and a vector plot of the mean in-plane motion of the flow field in different configura¬ 
tions are presented in the same figure. A wall-jet like boundary layer forms along the 
side walls of the curved pipes, which gradually strengthens with increasing the curva¬ 
ture. The secondary motion of the flow, which is merely 5% of the bulk velocity in the 
mildly curved pipe, reaches up to 15% in the strongly curved configuration at the same 
bulk Reynolds number. Due to the geometry-induced centrifugal force which acts always 
outwards, this side-wall boundary layer (SWBL) is subjected to a favourable pressure 
gradient as it moves towards the inner bend. Correspondingly, the fluid particles at the 
side walls accelerate to reach their maximum in-plane velocity magnitude close to the 
horizontal midplane (0 = 0). From this point, the in-plane flow decelerates approaching 
the stagnation point of the Dean vortices in the inner bend to lift up towards the outer 
bend. This behaviour can be clearly seen in the contours of the mean streamfunction 
{i.e. streamlines displayed in hgure[S]). 

The map of the mean secondary motion in the centre region of the pipe section largely 
differs for the various curvature configurations. In the weakly curved pipe the in-plane 
flow is accelerated towards the pipe centreline after reaching the stagnation point at 
the inner bend, whereas the mean in-plane motion accelerates back towards the side- 
walls in the strongly curved configuration. This generates a bulge region of minimum 
in-plane velocity magnitude just below the pipe centre of the case with n = 0.1. The 
iso-contours of the stream function in figure [5] display the resulting antisymmetric mean 
Dean vortices. The difference among the two curvature configurations is striking and 


ig tne fv 
/.I l201.3h . 


non-trivial to explain (INoorani et al\ 

The previous characterisation of the flow in a bent pipe was concerned with mean-flow 
quantities and as such could be also seen in laminar flow. To provide a clear view of 
the distribution of turbulent activity in the pipe section the normalised turbulent kinetic 
















Figure 4. Pseudocolours of axial velocity: (a) in a cross-section; (from left: k = 0, 0.01, and 
0.1 (b) in an azimuthal section for the straight pipe and the equatorial mid-plane for the mildly 

and strongly curved pipes. 


energy (TKE) defined as fc = (m'u')/ 2 , with the prime indicating the fluctuating velocity 
and implied summation over the dummy index i, is illustrated in figure [51 Compared to 
the straight pipe, the value of k in the inner side of the mildly curved pipe is smaller; 
however, the flow can still be considered turbulent for all azimuthal positions. On the 
other hand, the peak in the turbulent kinetic energy diminishes near the inner wall of 
the strongly curved pipe. The turbulence, hence, is significantly inhibited in this region 
and the flow is almost quiescent. It is clear from figure [6] (b) that near the pipe centreline 
the flow is more turbulent in the weakly curved configuration than in the straight pipe. 
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Figure 5. (Left-half) Pseudocolours of mean streamwise velocity of the carrier phase with 
iso-contours of the flow stream function, 'F. The black dot denotes the centres of the Dean cell. 
(Right-half) Pseudocolours of mean in-plane velocity of the flow (ue)^) with the vector 

plot of the in-plane motion, clearly showing the side-wall boundary layer and its near-separation 
from the wall; (a) the mildly curved pipe, and (b) the strongly curved pipe. 



Figure 6. Pseudocolours of turbulent kinetic energy (normalised with the respective (a) 

straight pipe; (b) curved pipe with k = 0.01; (c) curved pipe with k = 0.1. Note that, in order 
to offer a direct comparison no azimuthal averaging has been performed in the straight pipe 
case. 


Although in the majority of the pipe section with k = 0.1 turbulence is suppressed, it 
is increased near the region of the bulge visible in the secondary motion (c./. figures [5] b 
andO c). 


3.2. Particulate phase 

In the current study a total of seven particle populations with different SB are simulated. 
The density ratio is fixed to Pp/pf = 1000, and for each particle population, various SB 
are obtained by changing the radius of the particles. In order to imitate laboratory 
experiments the bulk Stokes number of each population is fixed for the various flow cases 
with different n. This corresponds to fixing the particle density and radius [i.e. fixing 
the physical particle) when varying the pipe curvature. The letter coding of particle 
populations and their corresponding parameters are presented in table [3] In the current 
paper, heavier or larger particles are referred to as Stp25, Stp50, and the smaller, 
lighter or less inertial populations are the Stpl and Stp5. The StplOO population is 
simulated mainly to investigate the asymptotic ballistic behaviour of the heavy inertial 
particles. 
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Case 

dp/Ra 

Sh 

’P{k=Q.Q) 

CT+ 

^^{k=0.0) 

‘^^(k=0.01) 

II 

O 

StpO 

n/a 

0.0 

n/a 

0.0 

0.0 

0.0 

Stpl 

3.727 X lO”'* 

0.0451 

0.0671 

1.0 

1.0503 

1.3020 

Stp5 

8.333 X 10"‘‘ 

0.2257 

0.1500 

5.0 

5.2516 

6.5095 

StplO 

1.179 X 10"® 

0.4514 

0.2121 

10.0 

10.503 

13.020 

Stp25 

1.863 X 10"® 

1.1284 

0.3354 

25.0 

26.258 

32.547 

Stp50 

2.635 X 10”® 

2.2569 

0.4744 

50.0 

52.516 

65.095 

StplOO 

3.727 X 10"® 

4.5139 

0.6708 

100.0 

105.03 

130.20 


Table 3. Parameters for the particle populations: Number of particles per population 
Np = 1.28 X 10®, and the density ration pp/pf is fixed to 1000 for all populations. The bulk 
Stokes number Stt, and particle diameter dp/Ra of each population is fixed while the curvature 
is varied in each flow configuration (k = 0.0, 0.01, and 0.1). 


3.2.1. Temporal convergence 

The particles are introduced into the velocity field of the fully developed turbulent 
flow with uniformly random distribution in radial r, azimuthal 9 and axial s directions. 
The initial velocity is set equal to the local flow velocity. To identify when the particle 
dispersion has reached a statistically stationary stage the instantaneous particle con¬ 
centration C (normalised with the mean particle concentration in the domain) in the 
viscous sublayer (y+ ^ 5; (1 — r/Ra) < 0.014) is monitored. Generally the near-wall con¬ 
centration is expected to increase due to the turbophoretic motion of the particles until 
it saturat es and th e proc ess reaches a statistically stationary stage. Recent analytical 
studies of Sikovskv ( 20141) suggests that a portion of these particles, which are drifting 
from the outer layer towards the near-wall region, may not possess enough momentum 
and engages into a slow and diffusional process before completely segregating near the 
wall. The time for this process is estimated to be inversely proportional to the particle 
Stokes number. Figure [7] shows the evolution of the observable versus the convective time 
of the simulation. As can be seen, in the straight pipe, after an initial increase, C levels 
off at t > 500i?a/wf) for large particles. For less inertial particles, however, this time is 
considerably longer as Stp5 reaches a steady state at about 1500i?a/ub while Stpl does 
not reach a statistically steady state within the simulation time. 

From figures [7] (b,c) quasi-periodic oscillations of the value of C can be observed in 
the curved pipes; these can be explained by the action of the secondary motion on the 
particles in the initial stage of the simulations. The initial random seeding leads to less 
particles near the walls and larger initial particle concentration close to the centreline of 
the pipes. Starting the simulation, the majority of these inertial particles is pushed by 
the Dean vortices towards the outer bend where they largely accumulate. The same Dean 
vortices transport the accumulated particles towards the inner bend via the side-walls 
boundary layers where they are eventually re-entrained back to the region of outward 
motion in the centre. The whole process is continuously modulated by the carrier phase 
turbulence until the particle concentration at the wall reaches an equilibrium state. The 
time period of the oscillations is longer for the lighter particles and shorter for the larger 
particles. This scaling suggests that the particles are mainly excited by the centrifugal 
force in this very initial stage. Nevertheless, it is clear that in these configurations the 
dispersed phase reaches the steady state at about t = 1500Ra/ub. Note that in the current 
plot C is only averaged spatially and no temporal averages were taken. Therefore, even 
after reaching a statistically stationary value, C still shows some fluctuations in time, 
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Figure 7. Instantaneous particle concentration C in the viscous sublayer {y^ < 5; 

{1 — r/Ra) < 0.014) normalised with the mean particle concentration in the domain as a function 
of time (convective time units): (a) n = 0.0 (b) k = 0.01 ,and (c) k, = 0.1. ■ StpO, • Stpl, 
-Stp5,. StplO, Stp25,-Stp50,-StplOO. 


due to the instationary near-wall turbulent velocity fluctuations. It is interesting to note 
that the level of these fluctuations is both dependent on the curvature k and the Stokes 
number. The largest fluctuations are observed for the heaviest particles in the strongest 
curvature, i.e. for those particles that are most ballistic, operating with the highest in¬ 
plane velocity. 


3.2.2. Instantaneous particle distribution 

The instantaneous particle distribution for various curvature configurations at the final 
time of the simulations t = 3000 is shown in figure [H The Stp50 population has been 
chosen for illustration purposes here. The projected front view shows that the particles 
are largely segregated at the wall in the straight pipe, as few particles remain in the outer 
layer. Some near-wall particle clusters can be seen in this figure, corre sponding to the pref¬ 
erent i al localisation of part icles with regards to near-wall dynamics (IMarchioli fc Soldati 
l2002t [^rdina et ^ 2012qI) . In comparison to the straight pipe, the amount of particles 
in the pipe core of the bent configurations is clearly larger. This is also true at the inner 
and the outer bend locations. At the same time, the core of the mean Dean vortices 
in the strongly curved configuration is essentially depleted of particles, and a notable 
plume of particles can be seen rising in the pipe centre towards the outer bend. Over¬ 
all, the particle configuration marks the Dean vortices that are normally not so clearly 
distinguishable in instantaneous snapshots of the flow (c./. figure |4] a). Furthermore, the 
intermittent (partially rc-laminarised) region in the inner bend is also characterised by 
harbouring only very few particles; it is clear that the particles follow the SWBL and the 
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Figure 8. (a) Projected instantaneous front views of particles (Stp50) distributed in various 
pipe cofigurations (from left: ft = 0, 0.01, and 0.1/ (b) A cut-away view of the particles (Stp50) 
distributed in the viscous sublayer for pipes with different curvatures. The direction of the flow 
is from left to right. 


region of highest TKE, as shown in the velocity and kinetic energy maps (c.f. figures [5] 
and [6]). 

Figure[8] (b) displays side views of the particle distribution in the viscous sublayer, again 
defined as all particles below ?/+ = 5, for each flow configuration. The Stp50 particles 
form clear particle streaks due to the preferential concentration of the particles close 
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to the pipe wall. These streaky structures are straight and randomly (in the azimuthal 
direction) distributed over the wall of the straight pipe. On the contrary, the elongated 
particle streaks in the bent configuration are curved and of helicoidal shape. This is 
clearly a consequence of the effect of the secondary motion. In fact, in straight channels 
and pipes, particles tend to preferentially align in elongated streaks correlated with slow 
departing fluid motions directed away from the wall and corres ponding slow streamwis e 
fluid velocity ; nam ely the ejection events of wall turbulence (jRouson fc EatonI 120011) . 
IPicano et al. (200^ showed that this tendency to accumulate in the ejection events is an 
essential condition for reaching a steady state particle concentration. It is, therefore, not 
surprising that the particle dynamics in the bent pipe shows the same phenomenology 
although being influenced and changed by the presence of the Dean vortices. 

In order to provide a clearer view of these near-wall particle streaks, the geometry 
of each pipe configuration is unwound azimuthally and unfolded. The Stp50 particle 
distribution in the viscous sublayer is shown in figure [S] in the s — 0 plane. For the 
straight pipe, k = 0, the particle streaky structures are homogeneously distributed over 
the pipe wall circumf erence. Also, in agreement with previous studies in channel flow 
(ISardina et al\v2Q\2(m a large-scale organisation is clearly visible which modulates the 
smaller-scale accumulation; these scales are of the order of the pipe radius, Ra- However, 
in the curved configurations the particle streaks are not only elongated in the axial 
direction of the pipe, but also inclined in the azimuthal direction to point towards the 
inner bend. In fact, these streaks are aggregated with a certain angle to resemble fish¬ 
bone like structures. The inclination of these helices appears to be three times larger 
in the strongly curved configuration compared to the mildly curved pipe. This can be 
explained by the fact that the in-plane motion of the carrier phase in pipes with /t = O.I 
is almost three times stronger than that of the mildly curved pipe. From this plane view 
(figure [9] b) it is clear now that the particles intensely accumulate in the inner side of 
the mildly curved pipe {9 ^ in/2) and form a large particle streak all the way though 
the pipe length. Interestingly the same position in the strongly curved configuration 
is totally void of these heavy particles as discussed above. A noteworthy increase in 
the particles segregation in the left half of the pipe with k = O.I is evident, which is 
changing to the right-half periodically in time at a comparably long time scale (not 
shown here). This changeover is reminiscent of the swirl-switching phenomena that is 
normally observed in spatially developing pip e bends and suggested t o be a result of up¬ 
stream effects or separation in the bends (see iKaloakli fc Orliilboi.lh . The fact that the 
particle concentration carries features of swirl-switching might indicate that these low 
frequencies are in fact inherent to strongly curved pipes even without upstream effects. 
A final observation from figure |9] is that the large-scale organisation seen for the straight 
pipe is absent for the bent pipes. The inclined particle streaks are much more uniform. 
This indicates that the duration a particle spends in the near-wall region is too short 
to be organised over spatial scales 0{Ra)- Similarly, for the pipe with strong curvature 
even the signature of the near-wall streaks is weaker, which could be associated with the 
stronger secondary flow. 


3.2.3. Particle trajectories 

Some typical trajectories of Stp5 and Stp50 particles within the last \2QQRa/ub con¬ 
vection time units of the simulation are shown for different curvatures in figure [TUI In the 
straight pipe these particles are already trapped in the viscous sublayer of the flow and 
remain there for this long period. Also, the particles hardly move in the azimuthal di¬ 
rections, and thus clearly sample one single low-speed region each. However, considering 
the same populations the trajectories are completely different in curved configurations. 
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Figure 9. An open-cut view of the particles (StpSO) distributed in the viscous sublayer for 
pipes with different curvatures; (a) straight pipe; (b) curved pipe with k = 0.01; (c) curved pipe 
with K = 0.1. The ordinates represent the azimuthal direction of the pipe and the abscissae are 
the pipe lengths with the flow being from left to right. 


In the case with k = 0.01 particles of Stp5 are erratically modulated by the flow tur¬ 
bulence unless they are trapped in the SWBL where they follow the secondary motion. 
On the other hand, the heavier particles of Stp50 are essentially following the mean 
Dean vortices performing a spiralling motion through the bent pipe. What is striking is 
the fact that both particles are spending a long time near the stagnation point of the 
Dean vortices in the inner bend once they are transferred by the SWBL and before being 
re-suspended towards the middle. In the case with k = 0.1 where the secondary motion is 
strengthened, the trajectories indicate that even the lighter particles are excited mostly 
by the act of Dean vortices. The trajectories of StpSO suggest that these heavy particles 
are hardly modulated by the flow turbulence for the majority of the pipe section. Once 
being trapped in the SWBL and descending towards the inner side, these particles are 
lifted up towards the middle section of the pipe in almost straight lines ignoring the car¬ 
rier phase in-plane streamlines; in particular the bulge feature discussed in figure [5jb). 
This event is followed by a series of wall collisions in the outer bend. These particle-wall 
collisions occur at an angle of incidence almost straight which might indicate a very large 
erosive effect (see lEdwards et al. 2001 ). 


3.2.4. Concentration statistics 

In order to obtain Eulerian statistics from the Lagrangian particle data in the straight 
pipe the whole computational domain is divided into wall-parallel axisymmetric slabs. 













































Figure 10. Some typical particle trajectories for (left) straight, (middle) mildely bent, and 
(right) strongly curved pipe, (top) The Stp5 and (bottom) the Stp50 population. The • indicates 
the beginning of the track and * shows the final time of the trajectory. Different colours indicate 
different particle tracks while the gray line sketches the circumference of the pipe section. All 
trajectories are 1200 convective units in length. 


The slabs in the viscous sublayer are distributed equidistantly with size of one particle 
diameter away from each other. Above this region ( 1 /+ > 5) the slabs are distributed 
non-uniformly according to the function: r = i?a(l ~ (sinh 777 )/(sinh 7 )) with 7 = 8 being 
the stretching factor. The vector rj consists of equally distributed grid points spanning 
the wall-normal direction. The resulting grid is more resolved near the wall than in 
the middle of the pipe. The curved pipe data is treated by an additional grid in the 
azimuthal direction. This mesh consists of cells (in the in-plane poloidal space) extended 
in the axial direction. A total of 2821 cells is used for the two-dimensional representation 
of the Eulerian statistics pertaining to the particle data in curved geometries: Ng = 91 
and Nr = 39, with 7 = 2; Ng and Nr are the number of cells in the azimuthal and 
wall-normal directions, respectively. Similarly, the statistical analysis at the equatorial 
mid-plane of the curved pipes is carried out with 100 wall-parallel slabs where 7 = 2 in 
the distribution function above and a constant cell width of 8 viscous units. 

The particle accumulation in various regions of the pipe section is quantified by means 
of the particle concentration C that is defined as the number of the particles per unit 
volume, normalised with the mean (bulk) concentration. Figure [TT] (a) displays this 
normalised mean particle concentration in logarithmic scale for each population in the 
straight pipe. While the Lagrangian tracers (StpO) are uniformly distributed (C = 1), 
the inertial heavy particles accumulate close to the wall and the maximum concentration 
appears at the wall due to the turbophoresis (see ICaporaloni et all 1975t Reeks! Il983l : 
Young fc Teeming 1997l l. Naturally, the particle radius serves as a cut-off of the profile 


in the vicinity of the wall. Despite averaging over a long period it is hard to say that these 
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Figure 11. (a) Wall-normal mean particle concentration profiles for various populations 

(normalised by the mean particle concentration in the pipe); ■ StpO, • Stpl,-Stp5, 

.StplO, Stp25, -Stp50,-StplOO. (b) The exponent of the power-law 

dependence C ~ (?/''')“ computed in the range of 1 < j/"*" < 5 (straight pipe). 


log-log concentration profiles show a linear behaviour. Though following [sikovs^ ( 2014 ). 
however, the exponent a of the power-law dependence C ~ (?/’'’)“ is computed in the 
range 1 < j/+ < 5 by a least-square fit. Figure [TT] (b) clearly shows that with increasing 
particle relaxation time this power-law exponent increases to reach a maximum at St~^ = 
25. This is of course consist ent with the the previous results of turbulent particle-laden 


chann el and pipe flows (see Young &: Teeming 1997 : Sardina et a/.l 12012^ IPicano et al 


20091) . which show that Stp25 particles are the most efficient in accumulating near the 


wall. Note that in the concentration plot, the lowest possible distance between particle 
centre and the wall is the particle radius, which increases with the Stokes number. No 
particle can be located below this minimum distance, which necessarily leads to zero 
concentration. 

Figure [12] displays the normalised mean particle concentration in logarithmic scale for 
mildly and strongly curved pipes for three different populations. Although particles in 
curved configurations also largely accumulate at the wall, the mean particle concentration 
is quite different from that of the straight configuration. As previously inferred from 
instantaneous snapshots (c./. figure ID a), the statistics clarify that heavy particles cluster 
densely outside the core of the mean Dean vortices, avoiding these vortical regions. As 
the strength of the Dean cells increases (in the case with k = 0.1) the concentration 
decreases as much as there are areas close to the vortex core where it is unlikely to find 
a particle at all. The cells in which no particles enter during a simulation interval of 
1500i?Q/Mb (the statistical time frame) are illustrated in white in figure [T2| (ej). The 
particle concentration in the inner bend stagnation point differs greatly for different 
curvatures. While the mildly curved configuration exhibits a high concentration in this 
region, in fact including the concentration maximum, whereas in the case with n = 0.1, C 
is very low (almost zero) in the same area. For both cases, this region grows in size with 
increasing the Stokes number of the particles. The strongly curved pipe also displays two 
high-concentration regions in the inner bend near the SWBL’s separation (c./. figures |5| 6 
andjOjc). These high concentration regions extend towards the pipe centreline and merge 
in the equatorial mid-plane to generate a particle confluence region. The bulge region 
visible in the velocity streamlines is not clearly represented in the concentration map 
of the heavier particles; however, smaller particles, e.g. Stp5, display an enlargement 
of concentration in this region. Both curved configurations present a relatively larger 
concentration near the outer side where the near-wall turbulent intensity of the carrier 
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Figure 12. Logarithmic representation {logio{C)) of mean normalised concentration map of 
the particles dispersion in bent pipes with k = 0.01 (left) and n — 0.1 (right). From top: Stp5, 
Stp25, and Stp50. The white cells indicate regions completely void of particles in the statistical 
time frame of the simulation. 


phase is relatively stronger. It is worth noting that the almost perfect left-right symmetry 
of these plots is clearly a sign of the statistical convergence of the present simulations. 

Figure [TS] shows the azimuthal distribution of particle concentration in the viscous 
sublayer of the curved pipes. In the mildly curved pipe, although the turbulence intensity 
of the carrier phase is higher in the outer half (7r/2 < 0 < tt), the near-wall concentration 
is higher in the inner half for almost all the populations except the smallest and largest 
particles (Stpl and StplOO). Closer to the inner bend, with increasing Stokes number the 
accumulation increases dramatically for Stp5 to Stp50 particles. The largest growth of 
the near-wall particle concentration from the outer bend towards the inner bend belongs 
to StplO population. 
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Figure 13. Azimuthal distribution of mean normalised particle concentration close to the wall 

< 5) of bent pipes with k = 0.01 (a) and k = 0.1 (b). + StpO,-Stp5, . StplO, 

Stp25,-Stp50, and- StplOO. The inner side and the outer bend of the curved 

pipes are located at 6/n = 3/2 and 0/tt = 1/2, respectively. 


The same analysis in the strongly curved configuration (figure [T3] b) shows that for 
heavier populations than Stp5, traversing from the outer bend towards the inner side, 
the concentration at the wall first slightly decreases and then increases again with the 
minimum around 0 « tt (see the inset of figure[T3] &). Approaching the inner bend farther 
the accumulation increases for smaller populations more dramatically. This behaviour 
continues up to a maximum in the particle concentration around 9 = 47r/3. Further 
down towards the inner bend the concentrations decrease sharply, clearly linked to the 
strong outward motion of the carrier phase and the attenuation of turbulence in this part. 
The rate of this drop in the concentration value has a direct relation with the inertia of 
the particles. 
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Figure 14. The ratio of concentration in the viscous sublayer of the bent pipes (k = 0.01 

-n; K = 0.1 ->!<) to that in the straight pipe, Cratio, for different populations. Black 

represents the ratio for the outer bend, light-gray for the inner bend and dark-gray shows Cratio 
at the horizontal plane. The close-up view is shown in the inset. 


The data points in figure [TT] shows the ratio of concentration in the viscous sublayer of 
the bent pipes (at horizontal and vertical cuts) to that in the straight pipe (Cratio) for 
various particle populations. The data for the bent pipes are obtained by averaging the 
particle data in near-wall sectors at horizontal and vertical cuts with (1 — r/Ra)'^ ^ 5 
and (RaO)'^ ~ 30.6. For both bent pipes, in the stagnation point of the Dean cells at the 
outer bend (black lines), Cratio of Stp5 particles is decreased by almost 50% compared 
to the Stpl population. Increasing particle inertia further, the Cratio continues to grow 
for both cases in the outer bend. However, the growth is much faster for the strongly 
curved pipe as Cratio ~ 1 for StplO and Cratio ~ 1-7 for StplOO particles while in the 
mildly bent pipe the Cratio of the StplOO population only reaches up to 0.8. The ratio 
of the concentration in the inner bend position, denoted by light gray, of mildly curved 
pipe to that in the straight pipe rises very rapidly with increasing particle inertia to 
reach a maximum for StplO particles where Cratio ~ 13 after which the value decreases 
again. On the contrary, the Cratio in this position for the strongly curved pipe is nearly 
zero for heavy populations. Notwithstanding, for Stpl particles the concentration in the 
inner bend for the strongly curved pipe is three times larger than that in the straight 
pipe. It is worthwhile to note that for the mildly bent pipe the inner bend concentration 
for Stp5, StplO and Stp25 particles is approximately 12 times, 32 times and 22 times 
larger that in the outer bend, respectively. Interestingly, in the near-wall horizontal cut 
(dark-gray lines) of the strongly curved pipe (where the magnitude of the secondary 
motion is maximum) the Cratio is almost constant at 0.6 for all the populations except 
Stpl. This is not the case for the mildly curved pipe as the trend of Cratio in the near¬ 
wall horizontal cut is identical to the outer bend profile, but with an almost 20% deficit. 

Further insight can be gained by measuring the wall-normal distribution of the particle 
concentration at the equatorial mid-plane of the curved pipes in polar coordinates (figure 
11511 . The behaviour in the strongly and mildly curved configurations is similar in the 
















Particle transport in turbulent curved pipe flow 


23 


middle of the pipe section (—0.7 < r/Ra < 0.7); namely the concentrations of the 
inertial particles are below the uniform distribution, but increase with larger particle 
inertia (except for lighter particles (Stp5-Stpl0) in the mildly curved pipe). Traversing 
towards the outer bend, the concentration of the particles generally reduces slightly. The 
accumulation of the particles in the inner bend is consistent with the observations in figure 
m However, an interesting feature of the particle concentration appears in the outer bend 
of the strongly curved pipe where there are series of peaks in the concentration map; 
these peaks are not a consequence of insufficient temporal averaging but rather a robust 
statistical feature of this flow. This can also be seen in figure[T2 (f). It is apparent that the 
magnitude of C and the radial position of these peaks is directly connected to the particle 
Stokes number. The reason for these concentration peaks is that particles hit the outer 
wall with considerable speed, and are thus bouncing from the wall (elastic collision). The 
incoming speed essentially determines how far the re-bounce will reach into the domain. 
If the incoming particles have a large-enough and quite uniform speed, such concentration 
peaks will be visible. The results show that the re-bounc e is stronger for heavier p articles. 
This phenomenon has also previously been observed bv iHuang &: DurbinI ( 20I2ll close to 
the outer bend of highly curved S-shaped channel laden with heavy particles and denoted 
as reflection layer. These authors explained the oscillations and the multiple peaks of the 
concentration in their study via a simplified model which shows that for large enough 
particle inertia the re-bounced ballistic individuals could have multiple velocities in the 
wall-normal direction. 


3.2.5. Force statistics 

In this section we try to elucidate the asymmetry in the concentration profiles in the 
equatorial midplane of bent pipes (observed in the previous figure) providing a compar¬ 
ison between two important forces contributing to the Stokes drag acting on particles 
in the vertical plane; namely, the centrifugal and turbophoretic force. It is worthwhile 
to note that the centrifugation of the mean Dean cells is absent in this plane due to 
symmetry. 

The distribution of the ratio of particle centrifugal acceleration Op^c = '^1^1 R to the 
magnitude of the total particle acceleration Op^tot = + “Sp + ®sp) averaged in time 

and streamwise direction in the equatorial midplane of both bent pipes is shown in the 
figure[T6l Here, Vs^ indicates the streamwise particle velocity while Or ^, and Osp denote 
particle acceleration in radial, azimuthal and streamwise direction of the pipe section. 
In the mildly bent pipe, the importance of the centrifugal force is increasing directly 
with particle inertia. Nevertheless, the maximum value for all the populations is in the 
outer-bend half of the pipe section between r/Ra ~ 0.5 and 0.7. These extrema traverse 
farther towards the outer bend with Stf,. For the heaviest population in the maximum 
point of the profile, Op^c is almost 90% of the total acceleration although {op^do-p,tot) 
remains smaller than unity for all the positions in the vertical cut of the mildly bent 
pipe. 

The importance of centrifugal acceleration is much more obvious in the strongly bent 
pipe (figure [16] (h)) as the ratio of {ap^dap,tot) is almost at unity in the majority of the 
pipe section. This value slightly reduces near the bulge (—0.6 < r/Ra < —0.2) due to 
the fanning effect of turbulence in this region. The large peaks at the outer bend are 
purely due to reflection laye r s obs erved in the concentration plots of heavy particles. 
As shown by iNoorani et ai ( 2015 ). these specular reflections change the direction of 
the radial velocity at the outer bend and expel particles from the near-w all region. The 


consec utive peaks of radial velocity near the outer bend (figure 6 (c) in iNoorani et a,I 


( 2015h 'l show how the radial velocity cancels due to the superposition of reflected particles 
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a 



Figure 15. Distribution of the mean normalised particle concentration along the equatorial mid¬ 
plane (vertical cut) of the curved pipe with k = 0.01 (a) and n = 0.1 (b). + StpO,. Stp5, 

-Stp25,-Stp50, and-StplOO. The inner bend is located at r/Ra — —1 and the 

outer side of the curved pipe is at r/Ra — 1. 


with negative radial velocity. Consequently, in these positions the magnitude of the in¬ 
plane velocities is drastically reduced and the only effective component is the streamwise 
particle velocity which will be directly translated into the centrifugal acceleration. The 
peaks in the vicinity of the outer bend wall in figure [16] (b) are directly corresponding to 
the troughs at the outer bend of radial velocity prohle shown bv iNoorani et al\ ( 20151) . It 
is also important to note that {ap^do,p,tot) is larger than one in 0 < r/Ra < 0.6 in hgure 
[T6l (b). This is only possible if a negative acceleration exists in this vertical cut to balance 
the particle centrifugation. Such decelerations are depicted also by Noorani et al. (2015) 
computing the radial acceleration in similar positions {c.f. figure 12 (d) therein). These 
negative radial accelerations are directly related to the lag of the particle centrifugal force 
behind the fluid cen trifugal force which st ems from different streamwise velocities of the 
particles and fluids. Noorani et al\ ( 20151) further show that the streamwise velocity of 
the particle phase is smaller than that of the fluid {c.f. figure 6 (a) therein). Therefore 
it takes some time for the inertial particles to reach the same local streamwise velocity 
as the fluid Us, leading to the phase lag and strong deceleration in R. 

Assuming Tp is small, the contribution of various forces to the in-plane particle dynam¬ 
ics in turbulent bent pipes is derived from the particle equation of motion in Appendix 
El Each term in the equations (lA 91) and (lA 101) can be computed from the carrier-phase 
statistics. Among the contributions to the turbophoretic drift of the particles (the sec- 
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Figure 16. Wall-normal distribution of the ratio of particle centrifugal acceleration ap,c to total 
particle acceleration magnitude ap,tot averaged in time and streamwise direction in the equatorial 

midplane of (a) the mildely bent pipe and (b) the strongly curved pipe. -I- Stpl, . Stp5, 

-StplO, Stp25, -Stp50, and-StplOO. The value of 1 is indicated by ■. 

The inner bend is located at r/Ra = — 1 and the outer side of the curved pipe is at r/Ra = 1. 


Olid and third terms on the right-hand side of equations (IA9|) and (lA 101) ) the only 
non-negligible term in the equatorial midplane is 

_ d{u’j.u'fl} 

^turb — 7^ • 

or 

Figure [T3 shows the absolute value of the ratio of this acceleration term to the mean 
centrifugal acceleration (oc = (ul/R)) denoted as Oratio- The ratio of the absolute value 
of the secondary-motion contribution to in-plane acceleration [asecondary) to Oc in the 
plane of symmetry of the pipes is also shown in the same figure (bottom). In the equatorial 
midplane, Osecondary IS the summation of the 6th to 9th terms of the right-hand-side of 
the equation (IA9I) in the Appendix. 

For both cases aturb is only significant in the vicinity of the walls. However, aratio in 
the inner and outer bend position of the mildly curved pipe is much larger than that of 
the strongly curved pipe. This is mainly due the reduction of turbulence activity in the 
highly curved configuration. For the strongly curved pipe, Oc is larger than aturb even at 
the outer bend where the flow is fully turbulent. The peaks of aratio near the inner bend 
of this case {k = 0.1) are related to stronger turbulence in the bulge region. Moving to 
figures [T3 (c) and (d), the secondary motions have only a negligible impact with respect 
to the centrifugal acceleration for the high curvature case (figure [17] d). For the mildly 
bent pipe, however, the effect is of the order of 10% compared to centrifugal acceleration, 
reaching a peak at 40% in the inner part of the bent pipe. 
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Figure 17. (Top) Distribution of the ratio of the main component of turbophoretic acceleration 
to the mean centrifugal acceleration aratio along the equatorial mid-plane (vertical cut) of the 
curved pipe with k = 0.01 (a) and k = 0.1 (b). The absolute value is plotted here. (Bottom) The 
ratio of the absolute value of the mean secondary-motion acceleration to the mean centrifugal 
acceleration in the plane of symmetry of the curved pipe with k = 0.01 (c) and n = 0.1 (d). 
The inner bend is located at r/Ra = — 1 and the outer side of the curved pipe is at r/Ra = 1. 


Comparing the two previous figures, it appears that in the outer layer of the mildly 
bent pipe in the vertical cut the centrifugal force is the dominant force, and this domi- 
nancy increases with particle inertia. However, at the walls the turbophoresis is still the 
dominant cause of the particle in-plane motion. In the strongly curved pipe, on the other 
hand, the centrifugal force dominates the particle transport even close to the walls. 


4. Concluding remarks 

In this paper, a series of new direct numerical simulations (DNSs) of the turbulent 
flow in straight and bent pipes have been discussed. The curvature parameter k varies 
from 0, 0.01 to 0.1. Together with the velocity and pressure fields, populations of iner¬ 
tial point-particles were advected. All simulations were performed at moderate Reynolds 
number, corresponding to Rbt = 360 in the straight pipe configuration; this Rcr is large 
enough for the flow to be not completely dominated by low-Re-effects, but yet the study 
is computationally feasible. A total of 120 million grid points was used for each simu¬ 
lation, and 0.9 million Lagrangian particles were advected in each case. The employed 
numerical method is a high-order spectral-element method coupled with spectrally accu¬ 
rate interpolation employed for the particle solver. The implementation and set-up were 
validated against existing channel-flow data. 










27 


Particle transport in turbulent curved pipe flow 


The DNS results show that even a slight curvature of the pipe geometry drastically 
changes the map of particle concentration due the cross-stream motion. The inner bend 
of the mildly curved pipe exhibits larger concentration than the outer side for a range of 
particle Stokes number. Enhancing the curvature alters the dynamics of the particles even 
further as the change in the concentration map is more pronounced. While the core of the 
secondary-flow Dean vortices and their stagnation points in the inner bend are almost 
void of any particles, considerable accumulations appear in the outer bend. Generally the 
heavier (i.e. larger) particles do not exactly follow the in-plane flow streamlines leading 
to a large collision rate at the wall at the outer bend. This generates a reflection layer 
with large concentration adjacent to the outer bend. This finding is consistent with the 
study of part iculate dispersion i n S-s haped turbulent channel at very strong curvature 
conducted by Huang fc DurbinI ( 2012l l. 

In addition to concentration map, the instantaneous flow structures were examined. 
It turns out that the particle streaks, commonly observed in canonical wall-bounded 
particle-laden flows, exist also in the near-wall region of the present curved geometries. 
However, due to the general tendency of the particles to follow the secondary flow, these 
particle streaks are inclined from the outer towards the inner bend, and as a consequence, 
their length and strength is weaker. In general, due to the presence of secondary flow, the 
particles are more active, and cover larger portions of the cross-stream plane. Particles 
are not only found in near-wall regions as in the canonical case, but pass through the 
pipe centre region. However, as discussed above, if the strength of the secondary flow 
becomes large, certain regions of the flow are depleted of particles, whereas others are 
oversampled. 

Another important feature is that the particle dispersion reaches a statistically station¬ 
ary state inside the curved configuration and there is no dependency of the inlet/outlet 
condition. The present case thus allows us to study particle dispersion in the presence of 
secondary flow without temporal transients. It is worth noting that the Dean vortices are 
locked in their position unlike the situation appearing e.g. in a curved channel or curved 
boundary layer where vortical Dean cells or Gortler vortices appear. In those cases, the 
spanwise homogeneity leads to uniform particles concentrations in that direction. The 
present case of bent pipe, as a first of its kind, provides a suitable configuration to study 
the effect of geometry induced centrifugal forces, the Dean cell centrifuges and the tur- 
bophoretic motion of the particles on the dispersion map in the pipe cross-section. The 
current study may provide a reliable base to conduct further theoretical studies on the ef¬ 
fect of the curvature on turbulent wall-bounded particle-laden flows. The data can also be 
employed as a bench-mark to facilitate modelling advancement for dispersed multiphase 
turbulent flow in wall-bounded complex geometries. 
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Appendix A. 

This appendix contains the derivation of different force contributions to the in-plane 
particle motion in turbulent bent pipes. In order to distinguish between the main three 
effects determining the particle motion, namely turbophoresis, preferential concentration, 
and centrifugal effects of the Dean vortices we start from the particle-dynamics equation 
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in the approximation of small relaxation time Tp. In this case the particle acceleration 
can be approximated by the fluid ac celeration dvn/dt D u/Dt where Vp is the particle 
velocity, and u is the fluid velocity ( Nowbahar et al. 2013ll : 


— Ud Tr, 


Du 


(Al) 


Here, Up is the fluid velocity at particle positions and D/Dt = d/dt + u-V is the material 
derivative following the fluid. 

Let us first briefly recall the implications of the previous equation in turbulent plane 
channel and straight pipe flows. In the case of channel flow, we indicate the wall-normal 
direction with y and apply averages to equation dAH) in time and in the streamwise 
and spanwise directions, denoted by an overbar. Projecting this mean equation in the 
non-homogeneous wall-normal direction we obtain 


(K. 

' dy ’ 


(A 2) 


applying Reynolds decomposition, / = f + f, where / is a ge neric scalar or vectorial 
field and the prime indicates the fluctuations around the mean ( Pood 2000ll gives 


'^Vp 


^Vp 


- T, 


du ' yU'y 

dy 


(A3) 


where the mean particle velocity is due to two different contributions: Uy^, that represents 
the preferential concentration, sampling of specific fluid events, and Tpdu'yU'y/dy the 
turbophoretic drift. For channel flow at steady state, the mean particle velocity on the 
left-hand side is zero, which implies a net balance between the preferential concentration 
and the turbophoretic drift, namely, = Tpdu'yU'y/dy. Since the sign of the derivative 
is negative, the the partic les tend to preferentia lly accumulate in the low-speed regions 
as shown among others in Sardina et al\ ( 2012afl . 

For the straight pipe flow, proceeding as for the case of the channel, i.e. averaging in 
time, streamwise and azimuthal {9) direction and projecting the resulting equation (lA 1[) 
in the radial direction r we obtain after Reynolds decomposition 




du''i 


U'U' — UnU 


■e^9 


dr 


(A 4) 


where the last term is negligible in the buffer layer (as shown in Picano et aL ( 2009h l 


Urp represents preferential concentration and Tpdu(,M5./dr the turbophoretic drift. The 
equation has the same structure as in channel flow, as the last term is almost zero, 
therefore the same conclusions as for the channel can be drawn. In particular, the particles 
tend to oversample the negative flow events once the statistical steady state has been 
reached. 

For the ben t pipe, we can m ake use of the orthogonal thoroidal coordinate system as 
introduced in German^ ( I982fl . We average only in time and streamwise direction (s) 
denoting this average with (•) and project the equation in the radial (r) and spanwise 
coordinate (9), obtaining 
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where k is the curvature of the bent pipe and the metric term hs is defined as 1 + Kr sin(d). 
The previous equations can be rewritten as 


.... dipij-Urffl 2/zg 1 , . 1 d{uj~UQ^ 

{VrJ = {UrJ - Tp -- Tp - j——{UrUr) -Tp-- 


P dr 
kcos{9) 


h^r 


hjc 


Ksin(d) 

(UrUe) + Tp - - - lUsUs) + Tp 


r 80 
(ugue) 


(A 7) 
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^ dr 
kcos{9) 


(ugug) + 


^ hsr 

kcos{9) 
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r 89 

(UrUg) 


^ h, hs ' ' ‘ r 

Applying Reynolds decomposition and rearranging we finally obtain 


{Vr^) = {UrJ - Tj 
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Here we assume the effect of the secondary vortex is only due to the mean flow. On the 
right-hand side of equations (lA 9|) and (jA 101) . the first terms are considered as preferential 
concentration, and the second and third terms as the contributions to turbophoretic drift. 
The forth and fifth term in equation (lA 91) and the forth term in equation (lA 101) are also 
considered as the contribution to centrifugal forces. Except for the last two terms in 
those equations which are negligibly small, the rest of the terms on the right-hand sides 
is directly related to the secondary-motion contributions to the in-plane particle motion. 
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